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Abstract 



Semileptonic decays of D-mesons into scalar hadronic states are investi- 
gated. Two extreme cases are considered: a) the meson decays directly into 
an uncorrelated scalar state of two two mesons and b) the decay proceeds via 
resonance formation. QCD sum rules including instanton contributions are 
used to calculate total and differential decay rates under the two assumptions. 
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I. INTRODUCTION 



Low lying scalar mesons are an old problem in hadron physics, see the review by Spanier 
and Tornqvist on scalar mesons in [I] and the literature quoted there. In a recent analysis 
Pll of D- and Z^-meson nonleptonic decays distinct signals for strong enhancements in 
S'-wave irir and Kir channels have been observed, reviving the interest for these states. In 
the following we shall refer to the signal in the tttt channel as a and in the Kir chanel as k. 
The enhancements can be well described by a Breit-Wigner-type resonance form, with the 
correct threshold behaviour 

PX ' BW{S) ~ (s-m 2 x ) 2 + m 2 x T x (s) 2 ' (1) 

where the subscript X stands for a or k . The correct threshold behaviour is guaranteed 
through the s-dependent width 



A 1 / 2 (s, ml, ml) m 2 
X 1 / 2 (m 2 X} m 2 a} ml 



ix\s) = Lax m/o/ 2 : ■ y z ) 



Here m a and m;, are the masses of the mesons in the decay channel and X(x,y,z) = x 2 + 
y 2 + z 2 — 2xy — 2xz — 2yz. In || the following parameters were found: 

1) for the S'-wave nir channel, m a = 777,5 = 777^, 

m a = 0.478 ± 0.024 ± 0.017 GeV; T 0a = 0.324 ± 0.042 ± 0.021 GeV. (3) 

2) for the S-wave Ktc channel, m a = mx, rrib = m^, 

m K = 0.797 ± 0.019 ± 0.042 GeV; T 0k = 0.410 ± 0.043 ± 0.085 GeV. (4) 

A theoretical analysis |3J of the non-leptonic decay D — ► an is in agrement with the 
experimental data. The semileptonic D-decays offer, in principle, much cleaner samples 
than the nonleptonic decays since there occur no problems connected with the presence of 
a third strongly interacting particle. In the non-strange channel important experimental 
information comes indeed from the analysis of the final state interaction of the pions from 
leptonic if-decays. For the case of strange mesons, the analysis of D decays can play a 
similarly important role. In a recent paper by the FOCUS collaboration || clear evidence 
was found that in the semileptonic decay D + — > K~n + n + v the D-meson does not decay 
exclusively into the hadronic vector channel, but that there is interference with a scalar 
contribution. In this paper we estimate the decay rates of semileptonic -D-decays into the 
scalar Ktt and tttt channels. We use the method of QCD sum rules which has been 
successfully applied to several semileptonic decay processes. For a recent review see [0] and 
the literature quoted there. 

Even if (broad) resonances in the hadronic decay channels exist, it is not clear whether 
they are due to interactions on the quark level, or if they are rather an effect of interactions 
in the purely hadronic channel, see for instance ||. Here a theoretical analysis can be very 
helpful. The semileptonic decay is supposed to occur on the quark level and therefore the 
decay rate should depend crucially on the direct coupling of the resonances to the quark 
currents. We consider two limiting cases 
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FIGURES 




FIG. 1. Schematic picture of the semileptonic vertex of a D-meson decay, c, q±, q 2 are quark 
lines, jfi denotes the weak current. X might be either a (broad) resonance or represent an uncor- 
rected two-meson state of mesons a and b. 



• the observed two-meson final state couples to the corresponding quark-antiquark state 
only through the a or k resonance. 

• the quark antiquark state couples to uncorrelated final mesons in an S'-wave. 

Since the signals we investigate in this paper are very broad the finite width has to be taken 
into account. Otherwise our analysis is based on the same principles and assumptions as 
the used in the sum- rule analysis of other semileptonic -D-decays [[|. We therefore refer to 
this paper for details and clarifications. 



II. KINEMATICS 



We investigate the semileptonic decays 

D -> Xlv t , (5) 

where X might be the a, k or an uncorrelated ttti or Kn pair in the S-state, as depicted in 
Fig.0. 

The semileptonic decay of a D-meson with momentum into a scalar state with total 



momentum px and invariant mass y/s = \jp\ is described by the two form factors of the 
matrix element of the weak current 

(X\j^\D) = (p D +Px)J+x(t) + (pD-px)„f-x{t) > ( 6 ) 

where t = (po ~Px) 2 - In the decay rate the form factor f_x is multiplied by the difference 
of the lepton masses and hence is negligible for both e and \x decays. 
The differential semileptonic decay rate is given by 

^T(M) = Gl\ V 1 2 3/2 , , pxQ 
ds dt 1927r 3 m|, D >J+x\> K , v > 

where Gp is the Fermi coupling constant and V cq2 the CKM transition element from the 
charmed quark to the quark q 2 . The total width is 
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I r{m D -VS? d 2 T(s,t) 

ds / dt , , . (8) 

(m a +m b ) 2 JO dS dt 

The spectral distribution in the invariant mass \/s of the hadronic final state is given by 

dT f(™ D ~^) 2 d 2 T(s,t) 

1 ' dt-jHr- (9) 



dy/s Jo ds dt 



III. SUM RULES 

The D meson in the initial state is interpolated by the pseudoscalar current 

j D (x) = c(x)z7 5 gi(x) , (10) 

where c is the field of the charmed quark and qi that of an up or down quark, summation over 
spinor and colour indices being understood but not indicated explicitly. The final hadronic 
state X is interpolated by the scalar current 

3x(x) = qi(x)q 2 (x) , (11) 

where q 2 represents a light quark field for X = a, and a strange quark field for X = k. The 
semileptonic decay rate is obtained from the time ordered product of the two interpolating 
fields in Eqs. fllQ|) and ([□]) and the weak current jjf = <?27^(1 — 7s )c 

TMiPht) =* 2 / d'xd'y(0\T[j D (x)j^(0) Jx (y)]\0)e^-^ . (12) 

In order to select the semileptonic decay rates into the lowest lying hadrons we insert 
intermediate states and obtain the following double dispersion relation in the phenomeno- 
logical side 

Ttx n {pWx,t) = iy ds D ds x (0\j D (0)\D)(D\j™(0)\X}(X\j x (0)\0) 

x tt^. k — h contributions of higher resonances . (13) 

sd -PdS x -tji x 



Introducing 



(0\ JD (0)\D) = A D p D (s D ); p 2 D = s D , (14) 
(X\ 3x (0)\0) = A xPx (s x ); Px = s x , (15) 



and using Eq. (|B|) we obtain 

rripnen/ 2 2 i\ 1 [ j , A A Pd{sd Px{s x 

T ax \PD>Px>t) = -o / ds D ds x A D A X 

71 J Sn — 



sd -Pd s x - p\ 
*{f+x(t, s D ,s x )(p D +p x )^ + f- x (t, s d ,s x )(Pd-Px)p) 
+ contributions of higher resonances . (16) 
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In the following we concentrate on the relevant form factor / + and introduce the factor 
T + which multiplies the vector (p D +px)^ 

T +x \PDiPxit) = -5 / ds D ds x A D A x 2 T J+x{t,s D ,s x ) 

n 2 J s D - p D s x ~ P x 

+ contributions of higher resonances . (17) 
For the D-meson the density Pd{$d) introduced in Eq. (|T4] ) is given by 

Pd{s d ) = ird(s D - m 2 D ) , (18) 

and we obtain 

A D = fDTn D /m c , (19) 

fo being the D-meson decay constant in the conventional notation. 
For the density px we use the two extreme ansatze mentioned above: 

• The quark-antiquark current j x in the J = + state couples to the a or K-resonance 
described by the Breit-Wigner distribution ([!]) 

Px(s) = p x ,Bw(s) . (20) 

• The quark-antiquark current jx couples to an uncorrected meson pair, the density 
being described by the density of two-particle phase space 



1oit z ° 



where m a , rrib are the masses of the mesons in the final state. 

The three-point function can be evaluated by perturbative QCD if the external momenta 
are in the deep Euclidean region 

^«(m c + mi) 2 , p 2 x < (mi + m 2 ) 2 , t < (m c + m 2 ) 2 . (22) 

In order to approach the not-so-deep-Euclidean region and to get more information on the 
nearest physical singularities, nonperturbative power corrections are added to the perturba- 
tive contribution 

*ST 04, A. «) = i\ l M S, x „„ «* C-A)£-A) + ? (W ( ° y> ' (23) 



The perturbative contribution is contained in the double spectral function a + x- The Wilson 
coefficients CV, multiplying the power corrections can be evaluated in perturbative QCD. 
The operators occur in the operator expansion of the time ordered product Eq. (p"2|); 
their vacuum expectation values, the condensates (0^), are introduced as phenomenological 
parameters. 
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In order to suppress the condensates of higher dimension and at the same time reduce 
the influence of higher resonances, the series in Eq. ( p3|) is Borel improved, leading to the 
mapping 

1 (-l) n p -™ 2 /M 2 

/<*W(JO. pr^-^W (24) 

Furthermore, we make the usual assumption that the contributions of higher resonances are 
well approximated by the perturbative expression 

— ds D ds x 7 ^tt , (25) 

with appropriate continuum thresholds %d and sqx- By equating the Borel transforms of the 
phenomenological expression in Eq. fllTD and that of the "theoretical expression" , Eq. (p3|) . 
we obtain the sum rule 

m r 7T J(m n .+m h ) 2 



f S »D fSQX 9 ,,,, 

7T JsD,th <>sx,th 

+ K +x (M 2 D ,M 2 x ,t) + K +I (M 2 D ,M 2 x ,t) , (26) 

where is some value below s x- in the zero width approximation we have of course 
sx = tti 2 x- K+x are the Borel transforms of the nonperturbative expressions due to the 
condensates and Kj is an approximation to the instanton contributions, which might be 
important in the scalar channel fT0 |. 

The decay constant fn and the coupling Ax defined in Eqs. fllil) , (|15D, and fll^) can also 
be determined by sum rules obtained from the appropriate two-point functions. Using the 
same procedure as described above we arrive at 

^(f D h ^(M 2 )) 2 e- m o/ M2 = f S0D dsa D {s)e- s ' M2 + K D (M 2 ) , (27) 

Tfl c ^ ' Jm^ 



and 



(Af eor (M 2 )Y(- r X dsp x {s)e- s ' M2 ) 

' ' ' V 7T J(m a +m b ) 2 ' 

S 0X ,., 2 



ds a x (s)e- s/M + K X (M 2 ) + Ki(M ) . (28) 
The analysis of the two-point function for the scalar mesons, and the explicit expressions 



for the functions occuring in Eqs. (pq) , fl27|) and fl28 ) are given in Appendices A and B 
respectively. 

The final sum rule for the form factor is obtained from Eq. ( p6|) by inserting for and 
A x the expressions f% eoT and A% eor of Eqs. (H and (||) 
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f+x(t,m 2 D ,s x 

in'n > . . 9. l r Sox 

7T J (m a +m b ) 

{ f$OD f$OX 

X 



fSOD fS0X 2 . , . 2 

/ ds D / rf^x S X , t)e~ s n/ M n e~ s ^ M ^ 

•> s D,th -> s X,th 

+ K +X (M 2 D , M x , t) + k +I (M*, M 2 X , t)) . (29) 



The radiative corrections for the scalar and pseudoscalar channels are known to be large 
ITJ . They are expected to be large in the three-point function too. By inserting the sum 
rule expressions for the two-point functions, Eqs. (|27|) and ( f28|) , in the denominator of the 
sum rule for the three-point function, Eq. fl29D , we expect, at least, a partial cancellation of 



these corrections [f[3,|13|] 



IV. EVALUATION OF THE SUM RULES AND RESULTS 



The sum rule Eq. ( E9h is evaluated in the same way as described in ||, and we only 
sketch the main steps of this evaluation. In the complete theory, the right hand side of 
Eq. ( p9|) should not depend on the Borel variables M 2 . However, in a truncated treatment 
there will always be some dependence left. Therefore, one has to work in a region where 
the approximations made are supposedly acceptable and where the result depends only 
moderately on the Borel variables . To decrease the dependence of the results on the Borel 
variables M 2 , we take them in the two-point functions at half the value of the corresponding 
variables in the three-point sum rules, i.e., in Eq.(^) we put 

M' D 2 = Ml/2 and M' x 2 = M 2 x /2 . (30) 

We furthermore choose 



M\ m 2 x 



Ml m 2 D — ml 



(31) 



We have checked that the results do not depend crucially on this particular choice. If 
the momentum transfer t to the lepton pair is larger than a critical value t cr , non-Landau 
singularities have to be taken into account |§ . Since anyhow we have to stay away from the 
physical region , i.e. we must have t -C (m c + mi) 2 , we limit our calculation to the region 
< t < t cr . In this range the t-dependence can be obtained from the sum rule (|29f) directly. 
It can be fitted by a monopole, and extrapolated to the full kinematical region. 

Since we do not take into account radiative corrections we choose the QCD parameters at 
a fixed renormalisaton scale of about 1 GeV 2 : the strange and charm mass m s = 0.16 GeV, 
m c = 1.3 GeV, the up and down quark masses are put to zero. We take for the non-strange 
quark condensate (qq) = — (0.24) 3 GeV 3 , for the strange quark condensate (ss) = 0.8 (qq), 
and for the mixed quark-gluon condensate (qcrGq)(qq) = ml(qq) with m\ = 0.8 GeV 2 . 

For the continuum threshold in the D-channel we take from || sqd = 6 GeV 2 . The 
standard value in the X channel would be s x ~ (% + 0.5 GeV) 2 , yielding s 0a ~ 1 GeV 2 
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and sok ~ 1.6 GeV 2 . As an additional condition we use the mass constraint from the two- 
point function, as described in Appendix A. 

We start with the decay D — > k£P^ and first consider the case where the scalar quark 
current ( |TTD couples directly to the ^-signal through a Breit-Wigner distribution (^). In 
Fig. 2 we show the different contributions to the form factor /+ K (0) in Eq. (p9|), as a function 
of the Borel variable M|>, using the continuum threshold Sq k = 1.6 GeV 2 . As the lower limit 
for the fiducial region in Mp we take that value of M|, where the perturbative contribution 
is one half of the total contribution. As upper limit we take the value = 15 GeV 2 , which 
is motivated in Appendix A. For such a high value of the Borel parameter the result is very 
stable but it is largely determined by the choice of the continuum model. The instanton 
contribution in the fiducial region is completely negligible; the five dimensional mixed con- 
densate is strongly suppressed compared to the three dimensional quark condensate. In the 
Borel variable of the V-channel the fiducial region corresponds approximately to the range 
1.7 GeV 2 < M 2 X < 5 GeV 2 . 

1.5 , . 1 . , . , . 1 . , . 1 




— 0.5 1 ' 1 1 1 1 1 ' 1 1 1 1 1 

3 5 7 9 11 13 15 

M D 2 (GeV 2 ) 

FIG. 2. Dependence of the form factor f+ K at t = on the Borel variable Mp. Here the decay 
is assumed to proceed through resonance formation and the state density is given by Eq. (p0[) . Solid 
curve: total contribution; long-dashed: perturbative; dashed: quark condensate; dot-dashed mixed 
condensate; dotted : instanton contribution. 

In the range < t < 0.5 GeV 2 no non-Landau singularities occur for our choices of the 
continuum thresholds. The momentum dependence of f +K can, in this i-range, be very well 
approximated by a monopole expression 

Mt) = , (32) 

and extrapolated to the full physical region. 
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For 1.4 GeV 2 < so K < 1.8 GeV 2 and for values of M|, discussed above, we find for the 
form factor at t — 

0.48 </ +K (0)< 0.55, (33) 

and for the pole mass, 

1.9 GeV <M P < 2.2 GeV. (34) 

The pole mass is considerably smaller than the mass of the charmed pseudovector-meson 
_D s i(2536) which would fit into the t-channel. In the D — > Kiv decay the pole mass also 
came out to be smaller |J than the mass of the strange vector meson D*(2114). 

In the limits of the Borel variables and the continuum thresholds discussed above we 
obtain for the total semileptonic decay width 

T(D -> Ktv) = (5.5 ± 1.0) x 10~ 15 GeV , (35) 

where we have used V cs = 0.97. The same calculation done in an - unjustified - zero width 
approximation would yield a total semileptonic decay width which is about 20 % larger. 

3 I . 1 . 1 . 1 . 1 



2.5 - 




0.6 0.8 1 1.2 1.4 



s" 2 (GeV) 

FIG. 3. Invariant mass distribution dT(y/s)/dy/s of the hadronic final state in the decay 
D — > (Kir)siv- Solid: resonance formation; dots: uncorrelated Ktt. 

The spectral distribution dT/dy/s, Eq. @, where y/s is the invariant mass of the tcK 
state, is given in Fig. 3, solid line. 

Next we investigate the same decay under the assumption that the scalar current in 
Eq. (|ll|) does not couple to a resonance, but to an uncorrelated nK pair in an S'-state, i.e., 
we use the ansatz (|2l|) for p K (s) entering the sum rules. In Fig. 4 we show the dependence 
on the Borel variable of the decay form factor /+«(0) for so K = 1.6 GeV 2 . Note that now 
the density (^1|) describes a two-particle phase space and therefore the dimension of f +K is 
different from the previous case. 



9 



15 

12 

7 — 9 
> 

CD 

a 

o" 

1 6 



3 





3 5 7 9 11 13 15 

M D 2 (GeV 2 ) 

FIG. 4. Dependence of the form factor f +K at t = on the Borel variable Mp. Here the decay 
is assumed to proceed through an uncorrelated TrET-meson pair in an 5-state, the density is given 
by Eq. (|2l|) . Solid curve total contribution; long-dashed perturbative; dashed quark condensate; 
dot-dashed mixed condensate; dotted instanton contribution. 

The masses of the pole fit to the t-dependence are practically the same as for the reso- 
nance case; for the form factor at t = we obtain 

5.0 GeV -1 /+ K (0) < 7.1 GeV" 1 . (36) 

The total width comes out for the case of an uncorrelated ttK pair in an S state as 

T(D -> {Kir) s lvt) = (3.7 ± 1.1) x 10~ 15 GeV. (37) 

The spectral distribution for this case is also shown in Fig. 3, with a dotted line. Although 
there is no resonance formation the resulting distribution shows a maximum at approxi- 
mately the mass of the k, which is an effect of the decrease of the total phase space near 
the kinematical limits. 

The evaluation of the decay D — ► alv follows exactly the same lines. Here the fiducial 
range in M|> is chosen according to the same criteria as before and goes approximately from 
M|) = 8 GeV 2 to 18 GeV 2 corresponding approximately to a range lGeV 2 < M 2 a < 2.3 GeV 2 . 
The continuum limit so<r was chosen between 1 and 1.6 GeV 2 . The resulting form factors 
and total decay width are, in the case of a resonance formation with a Breit-Wigner width 

0.42 < f +a (0) < 0.57 , (38) 

T(D -> <r£u t ) = (8.0 ± 2.5) x 10~ 16 GeV , (39) 
and for the case of two uncorrelated 7r-mesons in an 5-state 
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5 GeV" 1 < f +a (0) < 6 GeV" 1 (40) 

and 

T(D -> (7nr) s £P e ) = (4.5 ± 1.0) x 10~ 16 GeV. (41) 



3 



2.5 - 




FIG. 5. The invariant mass distribution of the hadronic final state of the rate in the decay 
D — > {ttk)s£v. Solid: resonance formation; dots: uncorrelated Kir. 

The spectral distributions for both ansatze are shown in Fig. 5. 

V. SUMMARY AND CONCLUSIONS 

The semileptonic decays of D-mesons into scalar hadrons are very similar to those into 
pseudoscalars. The same role played by the vector part in the pseudoscalar case is played 
by the weak pseudovector for the decay into scalars. Therefore, the theoretical expressions 
are very similar, only terms proportional to the mass of the strange quark change sign. This 
leads to a small reduction of the form factors compared to the decay into a pseudoscalar 
state. The main difference of the semileptonic decay rates into K and k is due to the 
different phase space. The decay into a mesons is suppressed by the small value of the weak 
(c, d) matrix element V c d ~ 0.225. The spectral distributions of the invariant masses of the 
Kit and ittt sytems are given in Figs. 3 and 5. Their maxima are well below the masses 
of the k and the a, and also the forms are quite different from Breit-Wigner distributions. 
The increase of the spectral ditributions is steeper and the fall-off substantially faster than 
for the corresponding Breit-Wigner forms. This is an effect of the total phase space in the 
semileptonic decay. 

If the scalar current does not couple directly to the resonance but only to an uncorre- 
lated meson pair, the decay rate is reduced by a factor two compared with the decay into 
resonances, but nevertheless an enhancement near 0.8 and 0.5 GeV is visible due to the total 
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final state phase space (see Figs. 3 and 5). With good statistics a discrimination between 
the two extreme cases should be possible. This would add valuable information on the na- 
ture of the intriguing low lying scalar resonances. In reality things might be complicated 
by a coupling of the interpolating current |Tl| on the quark level to a resonance as well as 
to a uncorrelated meson pair. In this case on has to construct a new density px and use it 
instead of the densities in Eqs. (|20| ) or (f21~f), in order to calculate the form factors fx(t). The 
procedure that follows is the same. 
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APPENDIX A: EVALUATION OF THE TWO-POINT SUM RULES 

In Fig. 6 we display the different contributions to A K , as a function of the Borel variable 
M' K 2 using the ansatz in Eq. (^) for p K (s) and s Qk = 1.4 GeV 2 . As expected the instanton 
contribution is particularly important for small values of the Borel variable. In spite of the 
weak dependence of the sum rule results on M' K 2 , it is well known that radiative corrections 
could be large in this case. Therefore, we do not attach great significance to the high stability 
of our results. 



0.1 



0.05 



> 

CD 

CD 



-0.05 

12 3 4 

M 2 (GeV 2 ) 

FIG. 6. Borel mass dependence of A 2 .. Solid curve: total contribution; long-dashed: perturba- 
tive; dashed: quark condensate; dot-dashed mixed condensate; dotted : instanton contribution. 
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In the zero width approximation one can obtain a sum rule for the mass by performing 
the logarithmic derivative of the right and left-hand sides of Eq. (p7| ) with respect to M~ 2 , 
assuming that fr> is independent of M~ 2 . In Fig. 7 we display the logarithmic derivative of 
the left and right-hand sides of the sum rule fl28|) , again assuming that A K is independent of 
the Borel variable. 

2 

1.6 

CM 

% 1.2 
O 
"o 

§ 0.8 



0.4 


12 3 4 

M 2 (GeV 2 ) 

FIG. 7. Borel mass dependence of the rhs ( dashed line) and the lhs ( solid line) of the sum 
rule generated by the logarithmic derivative of Eq. (|28| ) with respect to M~ 2 for the case of the k. 

For sq k = 1.4 GeV 2 there is a reasonable overlap between the two sides of the sum rule 
in the Borel window 0.9 < M 2 < 2 GeV 2 , and this defines the range for the Borel variable 
M'x in Eq. fl29|). Since all Borel variables are related by Eqs. (|30|) and (fTp), we obtain from 
the above range approximately 5 < M|, < 15 GeV 2 . 















\ 










\ 










\ 










\ 







































13 



0.06 



0.04 - 



o 



0.02 



- 



-0.02 



M 2 (GeV 2 ) 



FIG. 8. Borel mass dependence of A 2 . Solid curve: total contribution; long-dashed: perturba- 
tive; dashed: quark condensate; dot-dashed mixed condensate; dotted : instanton contribution. 



In the case of a, similar results are obtained and in Fig. 8 we show (for SQ a = 1.2 GeV 2 
and the ansatz in Eq. ( p0|) for p CT (s)) that the instanton contribution (dotted line) is even 
more important in this case, giving a very stable result as a function of the Borel mass. The 
quark condensate and mixed condensate contributions (dashed line) are zero since they are 
now proportional to the light quark mass, which we take equal to zero. 

In Fig. 9 we show the lhs and the rhs of the sum rule generated by the logarithmic 
derivative of Eq. (|28|) with respect to M~ 2 , as a function of the Borel mass, using the ansatz 
in Eq. fl20| ) for p a (s) and so a = 1.2 GeV 2 . We see that there is again a good overlap between 
the two sides of the sum rule in the Borel window 1.0 < M 2 < 2.5 GeV 2 . Since the mass of 
the particle is related with the square root of the lhs of this sum rule, from this figure we 
see that our result is compatible with m a = 0.5 GeV found in 0. 
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M 2 (GeV 2 ) 



FIG. 9. Borel mass dependence of the rhs ( dashed line) and the lhs ( solid line) of the sum 
rule generated by the logarithmic derivative of Eq. (p8|) with respect to M~ 2 for the case of the a. 



APPENDIX B: PERTURB ATIVE AND NONPERTURBATIVE 
CONTRIBUTIONS TO THE TWO- AND THREE-POINT FUNCTIONS 

In all this work we take into account the mass of the strange quark at most squared and 
neglect the mass of the light quarks. For the scalar meson X, we consider the particular 
case of k, since the a can be easily obtained from it by neglecting the strange quark mass. 
The perturbative contributions for the two-point functions denned in Sec. Ill are: 



m 



2\2 



od\s) 



(Bl) 



and 



(JkXs) 



2m 2 s ). 



(B2) 



The nonperturbative contributions including the quark and mixed condensates are 



K D {M 2 ) = -m c (qq)e- m ' /M2 



1 + 



m, 



2 ' m 2 . \ 



2M 2 



1 - 



2M 2 ) 



(B3) 



and 



K K (M 2 ) 



-m 2 JM 2 



(qq) + 



(88) 



+ 



m 



2M 2 



(qq) 



(B4) 



where we have defined (qg s a.Gq) = m^(qq). The instanton induced contribution is given by 
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n 



2m s m q 



M 2 z 



2 ..2 



X 



dX~ ____ga:-z 2 /4 

2 2/4 {x - z 2 /4) 2 



(B5) 



where z = Mp, and we have introduced the average instanton size, p, and the instanton 
number density, n, given by [[L(J 



n ~ -fm 4 , 
2 



p ~ -fm. 

F 3 



(B6) 



In Eq. ( p5|) m s and m q are the effective quark masses which we take to be m s = 400 MeV 



and m q = 300 MeV. 



The perturbative double spectral function is obtained by using the Cutkosky rules and, 
in the case of the meson k, is given by 



p +K (s D , s x ,t) 



-3 



6n 



! A 3/2 



where 



SXr 



-2s D s x t + m 2 c s x {sD - s x +t) + m 2 s D (sx — s D + t) + m c m s 
x (m 2 c (s D - s x -t)+ m 2 (s x -s D —t)- (s D - s x f + t(s D + s x 

X Q(s D - m 2 c )Q(s X ~ SXmin)Q(s Xm ax ~ S X ) , 



m 2 c (-m 2 c + m 2 + s D + t) + m 2 s D - s D t + (m 2 . - s D )yJ X(m 2 ,m 2 ,t) 



2mi 



(B7) 



(B8) 



and 



SXr 



m 2 (-m 2 + m 2 + s D + t) + m 2 s D - s D t - (m 2 c - s D )^j \(m 2 c ,m 2 ,t) 



2m? 



(B9) 



The nonperturbative contributions, including the quark and mixed condensates, which 
survive the double Borel transformation in p 2 D and p x are 



K +K (M 2 D ,M 2 x ,t) = (qq)e- m ^ M ° 



m c -m s 2 / 2 m c - m t 

m n m„ -. — 

2 ^ c 8M£ 

2m c — m s (4m 2 + m c m s — 2t)(m c — m s 

6M 2 D + 



2AM 2 D M 2 X 



m c — 2m s (m c m s — 2t)(m c — m s ) 

f2 ' 



6M 



D 



2AM 2 D M 2 X 



(BIO) 



and the instanton induced contribution is given by 

' ' dr ds du 



K +I (M 2 D ,M 2 x ,t) = 



tt z m s m q Jo 



(-fM 2 D - l)M 



x 



(1 + A 1S ) 2 ( 1S M 2 D - AM 2 X ) 2 1S < 



ru 



(s - r)(A - u) 



-p 2 (r+u) p -^m 2 p -a p 2 g sM^ 



(Bll) 



where 
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= s(yMl - 1) 
l + 4s 7 y ; 



S + (l + 4s 7 )t)Ml 



s + (1 + 4s 7 )t) - AM% 
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